A Simulation of Blood Cells in Branching Capillaries 
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The multi-cehular hydro dynamic interactions play a critical role in the phenomenology of blood 
flow in the microcirculation. A fast algorithm has been developed to simulate large numbers of 
cells modeled as elastic thin membranes. For red blood cells, which are the dominant component 
in blood, the membrane has strong resistance to surface dilatation but is flexible in bending. Our 
numerical method solves the boundary integral equationJ^'^ built upon Green's functions for Stokes 
flow in periodic domains.!^ This fluid dynamics video is an example of the capabilities of this model 
in handling complex geometries with a multitude of different cells. The capillary branch geometries 
have been modeled based upon observed capillary networks as shown in flgures [I(i!)] and [I(b)| The 
diameter of the branches varies between 10 — 20 fim. A constant mean pressure gradient drives the 
flow. For the purpose of this fluid dynamics video, the red blood cells are initiated as biconcave 
discs, as shown in flgure l(c)[ and white blood cells and platelets are initiated as spheres and 
ellipsoids respectively. The hematocrit as well as the average velocity vary in the different branches 
due to the nonlinear interactions of the cells. Physiologically, red blood cells enclose a concentrated 
solution of hemoglobin which is about flve times more viscous than plasma but for this simulation 
the viscosity of the fluid inside and outside of red blood cells have been taken to be the same. 
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Figure 1 : (a) Scanning electron micrograph of a microvascular corrosion cast of a capillary network 
within a terminal villus of the human placenta (with permission) P (b) Model capillary network 
used for this simulation. Cross sections vary between 10 and 20 /im. (c) Red blood cells are initiated 
as deformable biconcave discs. 




Figure 2: The flow of red and white cells and platelets (not in physiological proportions) in the 
model capillary branches. 
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